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We calculate the current- voltage characteristics and critical 
current j r J rray of an array of Josephson-coupled superconduct- 
ing wires. The array has two layers, each consisting of a set 
of parallel wires, arranged at right angles, such that an over- 
damped resistively-shunted junction forms wherever two wires 
cross. A uniform magnetic field equal to / flux quanta per 
plaquette is applied perpendicular to the layers. If / = p/q, 
where p and q are mutually prime integers, Ic" la, (f) is found 
to have sharp peaks when q is a small integer. To an excellent 
approximation, it is found in a square array of n 2 plaquettes, 
that /c" ay (/) ( n /?) 1//2 f° r sufficiently large n. This result 
is interpreted in terms of the commensurability between the 
array and the assumed q x q unit cell of the ground state 
vortex lattice. 



I. INTRODUCTION 

There has been much recent interest in Josephson junc- 
tion arrays with so-called long-range interactions. Such ar- 
rays typically consist of two layers of TV parallel supercon- 
ducting wires, arranged at right angles, in such a fashion that 
a Josephson junction is formed at each point where two wires 
cross (see Fig. Q). In this geometry, each wire in one layer 
is Josephson-coupled to all the wires in the other layer. If 
TV 3> 1, each wire has many "nearest-neighbor" wires to which 
it is coupled. Under these conditions, the thermodynamic 
properties of the array are accurateljj-idescribed by mean-field 
theory. Indeed, it has been showrJil that the usual mean- 
field theory for phase transitions in such arrays in a trans- 
verse magnetic fieldcl becomes exact in the limit TV — » oo, 
provided that each wire can be considered as a locus of con- 
stant superconducting phase. Experiments on finite arrays 
(TV « 17) have recently confirmed the accuracy of mean-field 
predictions.0 For very large arrays, the assumption that each 
wire is a locus of constant phase breaks down, and the mean- 
field theory has to be coHected to take account of the phase 
variation along the arrayfl This same effect causes the critical 
current of the array to saturate at a finite value which d«B|ends 
on the Josephson penetration depth, even as TV — > cxdBcI 

The exactness of mean-field theory is believed to have ah. 
analog in the dynamics of these arrays. Chandra et alEsu 
have studied a version of local dynamics equivalent, in the 
large-N limit, to the resistively-shunted junction (RSJ) model 
with Gaussian thermal noise and no applied current (they 
use a different numerical approach from that followed here). 
In the dynamical case, the effects of magnetic-field-induced 
frustration combine with those arising from the fact that each 
wire has many other "nearest-neighbor" wires with which it 
interacts. The result is a kind of glassy behavior, with many 
metastable states separated by large energy barriers, critical 
slowing down near the phase transition, and typically glassy 
aging effects. 



In this paper, we extend the RSJ model described above 
to include applied currents. Specifically, we calculate the 
current-voltage (IV) characteristics and the critical currents 
jjrray £ Qr arra y S Q f crosgec ] Josephson junctions, in typical ex- 
perimental geometries as shown in Fig. |IJ We assume that, 
even though TV S> 1, each wire can be taken as a locus of con- 
stant phase, and we study how / array is affected by a transverse 
magnetic field of magnitude / flux quanta per plaquette. Our 
principal result is that -f" ray (/) is extremely sensitive to /, in 
a manner related to the commensurability between the lattice 
and the vortex structure of the induced array supercurrents. 
A similar sensitivity is found in the array c-citical temperature 
T c (f ), as calculated by mean-field theory,Q and as confirmed 
in detail by experiments. It would be of great interest if simi- 
lar experiments could be carried out to confirm our calculated 

ir ay (/). 

We turn now to the body of the paper. In Section O we 
describe the geometry of our calculations and the model used 
to calculate the IV characteristics and 7" ray . Our numerical 
results are given in Section H, followed by a brief discussion 
in Section f VI. 



II. MODEL 

As shown in Fig. the array consists of two layers of par- 
allel superconducting wires arranged at right angles. A resis- 
tively shunted Josephson junction is assumed to form at each 
point where two wires cross. To determine J array , we assume 
that a constant bias current Ib is injected into one wire and 
extracted from a different wire. The IV characteristics are de- 
termined by calculating the time-averaged voltage drop across 
these two wires. We consider two configurations for current 
injection. In the first, the current is injected into a wire at 
the left edge of one layer and extracted from the wire at the 
right edge of the same layer (left-hand side of Fig. [I]) . In the 
other configuration, the input and output wires are centered 
in different layers of the array (cf. right-hand side of Fig. |l|). 
We also assume that a time-independent, uniform magnetic 
field B is applied perpendicular to the array, of magnitude 
B = /4>o/a 2 , where $o = hc/(2e) is the flux quantum, and 
a 2 is the area of each square plaquette. 

We assume that the phase of the superconducting order 
parameter is a isanstant on each wire, independent of position 
along the wire.u With this assumption, we can write down 
the equations determining the flow of current through the 
system. We assume an array of TV x M wires, and model 
the intersections between the crossed wires as overdamped, 
resistively-shunted Josephson junctions. The equations de- 
termining the time-variation of phases on each wire can then 
be written down using Kirchhoff 's Laws. For example, in the 
perpendicular current input configuration, the equations take 
the form 
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iElSjm = ^ [Ic Sm(V'j ~4>k - Akj) + hj + lL;kj] , 
k=l 

M 

I]35kn = 2_. {IcS'm(lpj -tj>k- A k j) + I k j + lL;kj] ■ (1) 
3=1 

Here Is is the constant bias current, which is injected into 
wire m in one layer and extracted from wire n in the other 
layer; (f>k is the phase of the fcth wire in one layer; and ijjj 
is that of the jth wire in the other layer. I c denotes the 
critical current of a junction formed at the intersection of 
any two wires (the junctions are assumed to be identical). 
The supercurrent flowing through the yth junction is then 
I c sin(tf)j ~<f>k~ Akj), where the quantity in parentheses is the 
gauge-invariant phase difference across the junction. We use 
a gauge such that the phase factor Akj arising from the vector 
potential of the magnetic field is Akj = 2-KkjBa? /$o-t9 Ikj 
denotes the normal current through the shunt resistance R in 
the junction at coordinates (ka,ja); it is given by 



h d 



- 4>k — A k j). 



(2) 



Finally, at finite temperatures, one must also include a 
Langevin noise current Ih;kj in parallel to the other currents. 
In the present calculations, we assume temperature T = 0, 
so that lL:kj vanishes.EJ The equations describing the par- 
allel current injection can be written down in an analogous 
fashion. 

Eqs. (|l|) can be conveniently expressed in a dimensionless 
form, such that current is measured in units of I c , and time 
in units of Ti/ (2eI c R), and then rewritten in a fashion which 
is more convenient for numerical solution. First, we rearrange 
Eqs. (^) to obtain 
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Of the N + M degrees of freedom in these equations, one may 
be chosen freely. We fix this degree of freedom by arbitrarily 
choosing 
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With this choice, and making use of the double-angle formulas 
for the sines and cosines, Eqs. (tBI) take the final form 
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The form of these equations is quite suggestive. Namely, the 
time-variation of phases in one layer depends only on certain 
averages over phases in the other layer. A similar result is 
known to hold for the averages of certain equilibrium functions 
of the phases in each layer, in the thermodynamic (i.e. large 
N and large M) limits. A consequence of this result in the 
equilibrium case is that mean-field theory is exactly valid for 
the thermodynamics of this system in the limit of large N and 
large M.tl 

Although we have not found an analytical solution to these 
dynamical equations, they are easily solved numerically for 
any N and M using a standard Runge-Kutta method. We 
have used a fixed-step-size, fourth-order Runge-Kutta algo- 
rithm, with time steps ranging between 0.01 and 0.1 dimen- 
sionless time units. We obtained the d.c. voltage V by aver- 
aging the instantaneous voltage over 100 time units. In order 
to obtain the IV characteristicsc-the input current was slowly 
swept from high to low currents.tj After each current step, the 
system was given 100 time units to reach steady state before 
the time-averaged voltage was computed. Zero-temperature 
calculations were carried out for two different N x N array 
sizes, N — 18 and N = 14, for each input configuration. 



III. RESULTS 

Figure [| shows the calculated IV curves for several fields 
in the parallel current configuration. The field is expressed 
in terms of the frustration f = Ba 2 /$o, the number of flux 
quanta per plaquette. By definition, the array critical current 
JJ rray is the largest bias current for which the time-averaged 
voltage (V)t = 0. By inspection of Fig. ^, J arlay is consider- 
ably suppressed by application of a finite magnetic field. If 
/ is a rational fraction denoted f = p/q where p and q are 
integers with no common factors, then the figure also shows 
that /c rray (/) is suppressed more at frustrations with larger 
values of q. 

These observations are made more quantitative in Figs. ^ 
and which show the calculated J" ray (/) for a wide range of 
fields. Figure ^ shows Ic" ay {f) for an 18 x 18 array in both 
the perpendicular (top) and parallel (bottom) input current 
configurations. In this case, /J rray (/) was obtained by di- 
rectly calculating the voltage across the array for a range of 
bias currents and recording the highest current with a zero 
mean voltage. In order to reduce the numerical noise, we 
have smoothed the critical current data as a function of frus- 
tration with a Savitzky-Golay filter £3 Clearly, Ic" ay (f) has 
sharp peaks near fractions / — p/q with small denominators 

q- 

Since it can be difficult to determine numerically the precise 
current at which (V}t — * 0, we have also calculated Ic" ay {f) 
by extrapolating the IV curve above J arlay , assuming the func- 
tional form 
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and determining m and j rarray by the least-squares fitting pro- 
cedure. The resulting data were also smoothed using the same 
Savitzky-Golay filter. This extrapolation algorithm, used to 
calculate the data in Fig. ^, considerably reduces the numer- 
ical noise near / array , and the smaller peaks in 7" ray (/) stand 
out more clearly than in the direct calculation. Using this 
procedure, we have calculated I array (/) for both an 18 x 18 
(top) and a 14 x 14 (bottom) array in the parallel current 
injection configuration. Clearly, there are peaks in IJ rlay (/) 
at all fractions / = p/q such that q < 16 in the larger array 
or q < 12 in the smaller array. We now discuss some possible 
reasons for this behavior. 



in Ref.H The proportionality between i array (/) and T c {f) is 
is the same as that which is believed to hold between criti- 
cal current and transition temperature in mare conventional 
two-dimensional Josephson junction arrays£j 

Finally, wa-|ppmment briefly on our results in light of the 
work of RefsEa. Their picture of frustrated arrays may well 
hold true even in the presence of an applied current. However, 
it is difficult to draw any quantitative conclusions from our 
own calculations about the energy barriers they discuss, other 
than those which are suggested by the field-dependence of 
the critical current. These results obviously suggest that the 
barriers are smallest for largest values of q, though the reasons 
for the analytical dependence (oc q^ 1 ^ 2 ) remain a matter for 
speculation. 



IV. DISCUSSION 

The present results show that / array (/) is a sensitive func- 
tion of the applied field /. At / = 0, the array is unfrus- 
trated and /J" ay (/) is maximized. At a rational frustration 
/ = p/q where p and q are integers with no common factor 
and q < N — 2, our numerical results suggest that height of 
the peaks in / array (/) approximately obeys the relation 

C'"'~/ /4 (s/P 2 ) 1/2 - (7) 

Here r is the number of flux lattice unit cells (assumed to be 
of size q x q) which can be inscribed in the array, s — rq 2 is 
the number of array plaquettes which are contained in all r 
of the q x q cells, and P 2 = (TV — l) 2 is the total number of 
plaquettes in the array (TV = M in a square array). This form 
can be made explicit by writing r = (intLP/g]) 2 , where int[x] 
is the greatest integer n < x. Then Eq. |t] can also be written 
as 

i a " ay ~ {q/P){int[P/q]f 2 . (8) 
In the limit of large TV, this result takes the even simpler form 

i a " ay ~ (n/q) 1/2 , (9) 

where n = TV — 1. The values predicted by eq. (?] are plotted 
as crosses in Fig. ||. Agreement with the numerical data is 
remarkably good. 

We can understand Eq. |^ as a measure of the commensu- 
rability between the field and the array. We expect that, just 
as in more conventional Josephson arrays with square plaque- 
ttes, the ground state energy configuration in a field f — p/q 
and no applied current has a flux lattice unit cell consisting 
of q x q plaquettes. A finite wire array at a given field can 
contain a certain number of such cells, with a number of pla- 
quettes left over. We may expect that the current required 
to depin such unit cells, and hence T" array (/), will increase as 
the number of flux lattice unit cells which can fit into the 
array becomes larger, and the number of "remainder" pla- 
quettes becomes smaller. Indeed our analytical expression, 
which fits our numerical data, expresses T" array (/) in this fash- 
ion. Note also that, according to Eq. for large values of TV, 
i" rly ~ q" 1 ^ 2 . Thus, in this regime, i alray (/) depends on q 
in the same way as the mean-field transition T c (f) reported 
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FIG. 1. Sketch of array geometry. The array consists of 
two layers of parallel wires, arranged at right angles. A resis- 
tively-shunted Josephson junction forms at each wire intersec- 
tion. There are two configurations for the current injection: 
parallel injection and extraction wires (shown on left); and 
perpendicular injection and extraction wires (on right). 

FIG. 2. Current-voltage curves at several frustrations for 
an 18 x 18 array with parallel input wire configuration. V de- 
notes time-averaged voltage between input and output wires. 



FIG. 3. The critical current Ic" ay (f) plotted versus frus- 
tration / for an 18 x 18 array in the perpendicular input 
wire configuration (upper plot) and the parallel configuration 
(lower plot). These data are obtained from a direct calcu- 
lation, smoothed using the Savitzky-Golay filterO The hash 
marks at top denote each rational frustration / = p/q with 
q < 16. Frustrations with p = 1 are denoted with bold hash 
marks. 



FIG. 4. i" ray (/) f° r an 18 x 18 array (upper plot) and a 
14 x 14 array (lower plot) in the parallel input wire config- 
uration. These data are calculated by a least-squares fit to 
the functional form, of Eq. |(| then smoothed using the Sav- 
itzky-Golay filter S3 Hash marks at top denote each rational 
frustration / = p/q with q < 16 (upper plot) or q < 12 (lower 
plot). Frustrations with p — 1 have bold hash marks. The x 
symbols denote the peak heights as predicted by Eq. ^, using 
the same, arbitrary normalization for each plot. 



4 




FIG. 1. Critical Currents of Josephson-Coupled Wire Arrays. J. Kent Harbaugh and D. Stroud. 




FIG. 2. Critical Currents of Josephson-Coupled Wire Arrays. J. Kent Harbaugh and D. Stroud. 
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FIG. 3. Critical Currents of Josephson-Coupled Wire Arrays. J. Kent Harbaugh and D. Stroud. 




FIG. 4. Critical Currents of Josephson-Coupled Wire Arrays. J. Kent Harbaugh and D. Stroud. 



